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Abstract 

The absolute and the relative quadrupole shape invariants and Kj, provide a model indepen- 
dent measure of triaxiality for /3-rigid nuclei. We will show that one can obtain (73 and from 
a small number of observables. The approximations which are made will be shown to hold within 
a few percent both in the rigid triaxial rotor model and the interacting boson model. The shape 
parameter is given for an exemplary set of nuclei and is translated into effective values of the 
geometrical deformation parameters [3 and 7. 

PACS numbers: PACS numbers: 21.10.Ky, 21.60.Ev, 21.60.Fw 
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I. INTRODUCTION 



One basic property of the nucleus is its geometric shape. Therefore, the nuclear shape, 
whether it is spherical, prolate, oblate, axially symmetric, or triaxial, is a key property of 
the ground state, as well as of excited states of the nucleus. Quantifying the nuclear shape, 
one usually turns to the well known geometric deformation parameters j3 and 7. These are 
deduced from a comparison of data with, e.g., the Davydov-Fillipov model of a rigid triaxial 
rotor [1] . This approach incorporates a major problem. A rigid rotor model cannot account 
for vibrations of the nuclear shape, which is a strong limitation. But, even if a model is able 
to describe also vibrations in the deformation parameters as, e.g., by the Bohr Hamiltonian 
[2], the geometric interpretation of the interacting boson model [3, 18], or the GCM [4], a 
second problem arises. In general the shape parameters 13 and 7 do not have fixed values, 
because the nuclei have in general not a rigid shape but they are vibrating. Thus it is useful 
to consider alternative parameters related to the shape of a nucleus, namely quadrupole 
shape invariants [5, 6, 7], which are model independent, and which are direct observables. 
In this paper we will discuss mainly the quadratic and cubic shape parameters ^2 and q^. 
We will focus on the relative cubic shape parameter — qzl^^"^ 1 which is independent of 
the nuclear radius Rq and the charge e. We will show that it is possible to obtain g2, <ii 
and K-i, with good accuracy from only few data. The cubic shape parameter K-i, is related 
to triaxiality and will be given for a variety of nuclei. Its connection to the geometrical 
deformation parameters will be discussed. 

Quadrupole shape invariants were introduced by Kumar [5] and widely used by Cline 
and co-workers, e.g. [6]. They are expectation values in a given nuclear eigenstate of higher 
order moments of the E2 transition operator, which is usually taken to be the quadrupole 
operator. Considering the ground state they are defined as 

92= e2(0+|(Q-Q)|0+) , (1) 

?3= e\^t\{QQQTM) , (2) 
54= e'{0t\(Q-Q)(Q-Q)\0t) , (3) 

where the dot denotes a scalar product and brackets denote tensorial coupling, Q is the 
quadrupole operator and e the elementary electric charge. Higher order moments can also 
be defined and are related to fluctuations in q^, but will not be discussed here. The moments 



2 



q2 and qz can be written in terms of averages of geometrical deformation parameters as 

q^ = e^Ql{P^) = e'QlPlff and ^3 = e^Ql{P' cos(37)) . (4) 

with 

Qo = ?,ZRI/{A'k) . (5) 

These quadrupole shape invariants can be renormahzed to the second order invariant 52 by 
17,8] 

^» = ^. W 
omitting the nuclear radius or the electric charge in this form. These quantities can in 
principle be obtained directly from data, but this is difficult in praxis because a large number 
of E2 matrix elements including signs is involved in expressions (1-3). This can be seen 
expanding the invariants g„ into sums over E2 matrix elements, which is shown here for q2 
and ^3: 

e'T.{^t \\Q\\n){n wQw^t) . (7) 

i 

^/l e'E(0^ II Q II 2+) (2+ II Q II 2+) (2+ II Q II 0+) . (8) 

An evaluation of q2 and q^ using extensive sets of experimental quadrupole matrix elements 
from multiple Coulomb excitation has been done for some nuclei by D. Cline and co-workers, 
e.g. in [6, 9, 10]. Of course, the existence of such extensive data sets is the favorable, but 
it is not the general case. Thus, there is great interest to obtain the shape invariants from 
more restricted sets of data. 



II. APPROXIMATIONS 



The basic idea is to invoke the Q-phonon scheme as has been discussed in [11]. This 
scheme was suggested by T. Otsuka [12], and was developed by a Koln- Tokyo collaboration, 
e.g. [13, 14, 15, 16]. The Q-phonon scheme implies that the wave functions of low-lying 
states are exhausted by only a few multiple Q-phonon configurations, where a Q-phonon 
itself is an excitation by the quadrupole operator. The 2i state in an even-even nucleus is 
dominantly a one-Q-phonon state. It was shown [14, 15] that the Q-phonon scheme holds 
with good accuracy for the lowest levels of collective nuclei. Here, as we will consider only 

3 



the very lowest states, we keep within the non-orthogonahzed Q-phonon scheme [16, 17], 
which will be shown to be sufficient for our purpose. The Q-phonon scheme gives a simple 
selection rule, namely, that an E2 transition between two states may change the number of 
Q-phonons in first order only by one, i.e. AQ = 1. Neglecting all Q-forbidden transition 
matrix elements with AQ > 2 gives the first order approximation. We will denote quantities 
given in this first order approximation by a superscript (1). This leads to a drastic truncation 
in the matrix elements needed in the expansions in Eqs. (7,8), e.g., q2 as given in Eq. (7) is 
approximated by 

q2 ^ g?^ = e2(2+||Q||0+)^ = 5(^2; 0+ ^ 2+) , (9) 

because transitions from the two-phonon 22 state or even higher-lying 2+ states to the 
ground state are Q-forbidden in first order. Eq. (9) reflects the well-known fact that in 
most even-even nuclei the largest part of the E2 excitation strength is concentrated in the 
flrst excited 2+ state. In the rigid rotor this B(E2) value is known to be directly proportional 
to the squared /3-deformation by Eq. (4). In the case of non-rigid /3-deformation Eq. (4) 
defines an effective deformation parameter (3^// or, making use of the approximation (9), an 
approximate /3gjj. 

Using the Q-phonon scheme in first order for one obtains 

Then, approximating the i^^a-parameter following its definition in Eq. (6) (n = 3) results in 

' {qiy/^ VlO (2+||g||0+) ' ^ ^ 

which is calculated from the ratio of the quadrupole moment of the 2^ state and its E2 
matrix element to the ground state. It turns out ,e.g., checking this approximation within 
the rigid triaxial rotor model or the lBM-1, that such a rude truncation of the sum given 
in Eq. (8) is not sufficient for a good approximation to K^, as we will show in sections III 
and IV. Therefore, we used a second order approximation, allowing in each term of the sum 
only one Q-forbidden matrix element with AQ — 2. Doing so, we derive a second order 
approximation for q^ as 



qP^qr'- = xl^e'[{2t\\Q\\0tr-{2t\\Q\\2t) + 2.m^ 

(12) 



which still includes only few E2 matrix elements. In the following we will always denote the 
second order approximation with the superscript appr. instead of (2), as it is the only one 
we use. Note that the approximation to q2 in second order approximation is the same as in 
first order, as a Q-forbiddcn matrix element would always appear squared and such terms 
are not included in this approximation, and we get 

qr"-- = and /3,7f ' = P^^^ . (13) 

Dividing q^^^' from Eq. (12) by q2^^', we get a second order approximation for that 
includes only four different E2 matrix elements, involving the lowest two excited 2'^ states 

and the ground state. 

A problem that appears is that the signs of the E2 matrix elements are needed, which 
are not known in most cases. Usually we know only the B(E2) values which are 

B{E2;J, ^ Jf) = e^JfWQUf . (14) 

This ambiguity in the signs can be avoided by using a relation between the signs of four 
matrix elements, suggested by Jolos and von Brentano [19]: 

sign((2+||g||2+)) = -sign((0+||g||2+)(2+||g||2+)(2+||g||0+)) . (15) 

This relation gives the relative phase of the two terms in Eq. (12). There is still an overall 
sign of K2, which is the sign of the quadrupole moment of the 2^ state, deciding between 
prolate and oblate deformation. Then, the second order approximation for K^^' is 



^r''- = J^sign(g(2+)) 



B{E2- 2+ ^ 2+) ^|B{E2■ 2+ ^ 0+) • B{E2- 2+ ^ 2+) 



\ B(E2; 2t Of) B(E2; 2t Of) 

(16) 

where we use an alternative but elegant definition of the squared quadrupole moment fol- 
lowing Eq. (14), 

1 35 
B{E2;2t - 2+) = - 6^(2+1 |g||2+)^ = — Qi2tr . (17) 

The approximation formula for given in Eq. (16) is the key result of this work. It 
allows to measure this observable directly and in a model independent way from only few 
data. These are four absolute B(E2) values, namely B{E2;2t Of), B{E2;2t Of), 
B{E2;2^ 2f), and B{E2;2f 2f), and the sign of the quadrupole moment of the 



2i state, which we consider as a fifth observable. This new method to determine is of 
particular importance because is closely connected to the triaxiality of nuclei, i.e. to 
7-deformation. For axial symmetry — —1 for prolate (7 = 0°) and — +1 for oblate 
(7 = 60°) nuclei, while A'3 drops to zero at a maximum triaxiality of 7 = 30°. This holds 
for geometrical models like the Davydov-Fillipov triaxial rotor model, as well as for the 
dynamical symmetries of the IBM. One major difi^erence between these two models is that 
the IBM describes non-rigid /3- and 7-deformation, e.g., in the U(5) vibrational limit and the 
0(6) hmit of 7-soft nuclei, in both of which vanishes. In the SU(3) and SU(3) dynamical 
symmetries of the IBM, which correspond to the prolate and oblate axially symmetric rigid 
rotors, respectively, the same values for are derived as in the geometrical model. In the 
following we will check to which extent 

K, ^ Kt^- (18) 
holds, using as a test the rigid triaxial rotor model of Davydov and FiUipov and the IBM-1. 



III. THE RIGID TRIAXIAL ROTOR MODEL 

The Hamiltonian of the Davydov-Fillipov rotor model is 

Hgeo = A^JI + A2JI + , (19) 

where Jn are the projections of the spin J on the three symmetry axes, and where the 
parameters are connected to the moments of inertia by 

n 



(20) 



26^ 

The moments of inertia can further be written in terms of the geometrical deformation 
parameters (5 and 7, 

e, = 4S/3^sin2(^7-f ^) . (21) 
The E2 transition operator is given by 

r(E2),eo = eg2M = eQo (3 cos(7) + V2 {D% + Dl^) sin(7)] , (22) 

where the are the Wigner-D- matrices and Qo is given by Eq. (5). We stress that this 
model with rigid (5- and 7-deformations is applicable to only a limited number of nuclei. 
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Nevertheless, apart from our discussion of the ground state deformation, the model is often 
applied to highly excited and strongly and also super-deformed bands as well, for which, in 
principle, our approach of X-parameters may also apply. 

In our calculations we vary the parameter 7 over the range (7 G [0°,30°]), covering the 
range of prolate axially symmetric and triaxial structures inherent to the model. The results 
for 7 e [30°, 60°] are fully symmetric to those given and thus omitted. The choice of (3 is 
arbitrary, as in the rigid case is independent of (5 and is given by 

/33cos(37) 

^3 = (^2)3/2 = - co«(37) • (23) 

In a similar way one defines an approximate deformation 7"^^^- from K^^^' by 

K^PP'- = -cos{2,-i''PP'-) . (24) 

In order to avoid the division by zero, we use the ratio 

as a measure of the quality of the approximation (18). The solid curve in the left panel of 
Figure 1 shows the quantity i?^^ versus the deformation parameter 7, calculated numerically 
using the code Davidov [20]. In the axially symmetric limit at 7 = 0° the approximation 
is exact. This also holds for the case of maximum triaxiality at 7 = 30°, while i?^^ is 
small for all intermediate cases with a deviation from one of 8% in maximum. The dashed 
curve represents the same calculation, but using the first order approximation ^3^^ from Eq. 
(11). The deviation of the first order approximation is clearly much larger with a maximum 
of about 30%, showing that the use of a second order approximation is unavoidable for 
transitional nuclei. On the right hand side of Figure 1 the corresponding absolute deviation 
of 7"P?"'- derived from Eq. (24) from the real 7-valucs in the model is shown as a solid curve. 
The maximum deviation is below 3.5° at 7 ~ 15°. Again, the deviation is much larger using 
only the first approximation, given as a dashed curve.. 

IV. THE INTERACTING BOSON MODEL 

Now, we check the quality of K^^^^' in the IBM-1, within the Extended Consistent Q 
Formalism (ECQF) [21, 22] and the Hamiltonian [8] 

HiBM = (1 - C) rid - ^ • (26) 



depending on only two structural parameters, ( and x, and omitting an overall energy scale. 
The E2 transition operator in the ECQF is chosen to be proportional to the quadrupole 
operator in the Hamiltonian, 

T{E2)iBM = e^g^ = eB[(s+d + d+s) + x{d^d)] , (27) 

with an effective boson charge cb, and rid — (d'^d) is the boson number operator. Varying 
the values of ( and x over the full range of symmetries {( e [0,1], x £ [— '\/7/2, -\/7/2]), 
one covers the dynamical symmetry hmits of the IBM, namely U(5) (C = 0, x), the prolate 
(oblate) SU(3) (SU(3)) (( = l,x = tV7/2), and 0(6) (( = l,x = 0), as well as the 
transitional structures in between. In analogy to Eq. (25), we define the ratio 

^ 1 + \K^"^'\ 

^™ " ~T+]k^ ' ^^^^ 

which has been calculated over the full parameter space using the code Phint [23] . Again, 
j^appr. -g (^jgf^^g^j \yy jT;q (15) The results are shown in the top part of Figure 2 for = 10 
bosons. They are given for x < 0, because the results for ±x are fully symmetric, as the 
change in sign is equivalent to the symmetry transformation d — —d (keeping s — > s). The 
use of positive x values corresponds to the choice of 7 > 30° in the geometrical model. 
Deviations of K^^' from the exact K3 values are small in all cases, the deviation of R^bm 
from 1 is below 7%. For comparison. Figure 3 shows R^b^m ^ which is defined analog to Eq. 
(28), but where the first order approximation, i.e. Eq. (11), is used. Like in the geometrical 
model it is seen, that the deviations from the exact value of K-^ are much larger in the first 
order approximation. Thus, in general it is necessary to use K'^^' from Eq. (16). 

The deviation RfsM peaks in a region around SU(3). A reason for this behavior is found 
by a close look at this region. The middle part of Figure 2 shows the values of and K^^' 
on the U(5)-SU(3) (left) and 0(6)-SU(3) (right) transition legs. It is obvious that the 
maximum deviation of R^bm from unity appears in those regions, in which changes most 
rapidly. These are exactly those regions that arc connected to the shape/phase transition 
between spherical and axially symmetric nuclei, or between prolate and oblate deformations, 
as discussed, e.g., in [11, 24, 25, 26, 27, 28]. This means that in the IBM the approximation 

j^appr. 

misses the exact value of somewhat when leaving the rotational limit. However, 
overall deviations of R^bm from unity are small and the approximation (18) is well fulfilled. 
From comparison with the geometrical model an effective 7-deformation can be defined 
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[8] from Ks by 

- ^"^^^,^3/2 = - cos(37e//) , (29) 

and an approximate value 'y'^fY' can be defined analog from K^^'. The differences between 
the exact and the approximate 7- values, 7e// — %fY included in the bottom part of 
Figure 2 and show good agreement. The deviation of %ff^' from 7e// is always smaller than 
2.5°. 

These effective 7- values are not and cannot be equivalent to those given by Eqs. (23,24), 
because is not generally independent of /3-deformation and fluctuations in f3 occur, es- 
pecially for vibrational nuclei. Moreover, in case of rigid (3 (on the SU(3)-0(6) transitional 
line) is a measure of (cos(37)), while in case of rigid 7 it is a measure of 
The effect of a /3-vibration is only effectively taken out in the translation to the geometric 
model by Eq. (29). However, if fluctuations in /3 are small, which is the case past the phase 
transition towards deformed nuclei (typically for ( > 0.6), a factorization of the averages 
over P and cos (37) should work, and we can assume 

(/3=^cos(37)) = (/3'^)(cos(37)) and = 1 , (30) 

making 7e// comparable to the geometrical 7-deformation. 



V. K3 FOR VARIOUS NUCLEI 
A. Direct measure of 

For the two considered models we have shown that K^^^^' is a good approximation to 
the value of the cubic shape parameter K3. Thus we assume this to hold also in other 
collective models such as the GCM or the Bohr Hamiltonian. Only few observables have to 
be obtained in order to derive K^^^^' , namely the hfetime of the 2^ state, the hfetime and 
the branching ratio of the 22 state, and the quadrupole moment of the 2i state. Besides the 
modulus of K3 also its sign is interesting, which is obtained from the sign of the quadrupole 
moment of the 2f state. This quadrupole moment itself is not easy to obtain, therefore it 
is a challenge to measure triaxiality. Especially for vibrational or 7-soft nuclei, where the 
quadrupole moment is small, high quality data is needed. Thus, an approximate value of 
Ks is so far known for a number of nuclei in or near the valley of stability only. For a set 
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of nuclei that belong to various symmetry regions the Xg^^'^ -parameter has been calculated 
from tabulated data. The results are given in Table I, together with effective 7-deformation 
parameters derived from Eq. (29), and effective /3-deformations from Eqs. (4,9). 

Typical rotational nuclei like the heavier Gd or Dy isotopes show values around -1 as 
it is expected for prolate deformed axially symmetric shapes. Also ^^^Sm and ^^^Gd, which 
are attributed [25, 29] to be close to the critical point symmetry X(5) proposed by F. lachello 
[24], show this value. The ir4-parameter obtained from 54 of Eq. (3) and Eq. (6), which 
can be approximated in a similar way [11, 19], gives a direct measure for /3-softness. One 
finds that Kf^- = 1 for /3-rigid nuclei and Xf^''- ~ 1.4 for vibrators. For ^^^Sm and ^^^Gd 
one finds K^^' values of 1.02(3) [30] and 1.088(26) [29], respectively, which, in combination 
with K^^' , meets the expectations for the vibrator to well-deformed rotor transition, more 
on the rotational side of the phase transition. Especially ^^^Sm seems to be on the rotor side 
(where K^^' — —1 and Kf^' — 1) of the phase transition, which seemingly conflicts with 
the interpretation of this nucleus as being close to the phase transitional point. This may 
be related to the systematical error made in the approximations for K^, which maximizes 
exactly in the transitional region in the IBM, which may be reflected also in data. However, 
the systematical error made in the determination of should be smaller in that region [19] 
and a problem remains. 

The values of the Os isotopes show an evolution from the axially symmetric rotor 
towards 0(6) symmetry with a maximum effective triaxiality of 7e// = 30°. Note, that here 
one talks of effective 7-deformation, as the nucleus does not have a rigid triaxiality. The 
more vibrational Pd and Cd isotopes show moderate values of with relatively large errors 
due to the quadrupole moments. Non-zero values are not a contradiction to a more U(5) 
like structure as they may emanate from finite N effects (see [8]). 

An surprising conflict appears for ^^^Pt, which is usually taken as a prime example of 
0(6) symmetry [31], as well as for the neighboring ^^^Pt. Both nuclei show rather large, 
positive quadrupole moments [9] and thus have quite large values of K^^' . This shows 
that they are on the side of oblate deformation, with a considerable deviation of K^^^^' from 
the expectation value, ^^3(0(6)) = 0. Other observables like the branching ratio of the 
22 state or energies agree much better with the 0(6) predictions. Therefore, values of 
derived from an IBM fit (see below) agree much better with K3 — 0. One cannot argue 
that this deviation is due to the dependence on /3-fluctuations (compare Eq. (29)). On the 
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SU(3)-0(6) transition line, no /3-fluctuations are allowed, and indeed, the shape invariant 
approximately equals 1 (see [11]) for both nuclei, which pinpoints /3-rigidity. Again, 
this may be related to the maximal systematical error close to 0(6) seen from Figure 2. 

But, even if the value of K3 is overpredicted from the approximation, a deviation from 0(6) 
remains. However, wc want to stress that these values, e.g. 'feff = 42° instead of 7e// = 30° 
for ^^^Pt, still indicate a strong triaxiality. It is only the quantitative value of jeff which is 
in doubt. 

We stress that if one uses only the first order approximation, the value of K3 is missed 
for transitional nuclei hke the Os isotopes, for which the transition 23" Of is sizeable, 
e.g., K^^^^' = —0.7(1) for ^^°0s in the first order approximation, underestimating triaxiality, 
while the second order approximation gives K'^^' = —0.35(9). 

Note, that for the electric quadrupole moment of the 2f state, Q{2i), is usually not easy 
to access experimentally. A new relation was found [11], however, which gives a new way to 
approximately determine Q{2f) or B(E2;2i — > 21"), respectively: 

B{E2; 2+ ^ 2+)"^'^"^- = B{E2; 4+ ^ 2+) - B{E2; 2+ ^ 2+) . (31) 

This relation may be used to obtain the quadrupole moment as an input for K3, but it 
will give a large uncertainty especially for vibrational or 7-soft nuclei, which have a small 
quadrupole moment. So far the relation (31) was only checked in the IBM [11]. Figure 4 
shows the deviation 



B{E2;2t ^2t)+B{E2;2t ^2t) ^ ' 



from the real value calculated within the rigid triaxial rotor model. Also in the geometrical 
model the agreement is good. 



B. Fit procedure for 

In cases where not all of the needed data are present, one may follow another procedure, 
fitting parameters of a model to the available data for one nucleus, and calculating from 
the model. Here we used the simple two parameter Hamiltonian of the IBM given in Eq. 
(26). The two parameters were fitted to the energy ratio 

i?4/2 = E{At)/E{2t) (33) 
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and the B(E2) ratio 

B{E2- 2+ ^ Qt)lB{E2- 2+ ^ 2+) , (34) 

that are sensitive to changes in structure over wide parameter regions. For the reproduction 
of the energy ratio an error of 2% was allowed, while for the B(E2) ratio the experimental 
errors were taken into account, resulting in an allowed parameter range of C, and x, in which 
takes various values within a certain range, with an upper a lower limit. In Table I we 
denote values of obtained from the fit as -ft'/**, and give the upper and lower limits allowed 
from the experimental errors. These values can be compared with the measured K^^^^' . The 
values agree reasonably well in most cases, considering the simplicity of the Hamiltonian 
and the arbitrary choice of the two observables used in the fit. Note, that other observables 
can be used for the fit. But in some cases the simple Hamiltonian used cannot describe all 
features of a given nucleus, as, e.g., for the Pt nuclei. Therefore, the Hamiltonian (26) may 
be extended or another model used. 

VI. CONCLUSIONS 

To conclude, we discussed measures of triaxiality. In this respect we considered in par- 
ticular the absolute and relative cubic shape parameters ^3 and K3. The approximative 

j^appr. 

was introduced as a direct, model independent observable, which is, if P is rigid, a 
measure of triaxiality, while K^^' is more general an observable in all structural limits and 
the regions between them. K^^' was shown to be a good approximation to the exact value 
of and can with good accuracy be obtained from only four matrix elements, or from 
four B(E2) values, one of them equivalent to the modulus of the quadrupole moment Q{2i), 
and the sign of Q{2]^). This, manifested in Eq. (16), is the key result of this work. This 
accuracy of the approximation was checked for the IBM and the rigid triaxial rotor model, 
and the need of the second order approximation within the Q-phonon scheme was shown in 
both models, especially in transitional regions. Effective values of /3- and 7-deformation in 
the ground state, derived from 02"^"^' and K^^^^', respectively, have been deduced from data. 
For vibrational nuclei geometrical deformation parameters cannot be given, while q2 and 
are always well-defined properties of the ground state. Finally, we proposed a way how to 
derive K^^^' from a model fitting data. 
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FIG. 1: R^^g calculated for all values of 7 (solid line on the left hand side), showing that the second 
order approximation Eq. (18) holds well in the rigid triaxial rotor model. The approximation does 

not seriously change the value of 7 (solid line on the right hand side). The dashed lines give the 
values derived from the use of only the first order approximation. 
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TABLE I: Approximate iir3^^'''-values, the effective approximate /9- and 7-deformation parameters 
derived from our approach are listed for a set of nuclei. For /3-deformations, errors are omitted as 
they are in the order of per mil or smaller, and the systematic error made by assuming Rq = 1.2 fm 
for the nuclear radius is presumably larger. The last two columns give upper and lower limits for 
the value of fitted to the observables (33,34) as described in section VB. 









appr. 


K 
upper 


fit 

T 
J 

lower 




-0.97(5) 


0.339 


4(4) 


-0.86 


-0.98 


is^Gd 


-0.95(6) 


0.349 


6(6) 


-0.83 


-1.00 


leoGd 


-0.96(3) 


0.351 


5(3) 

V / 


-0.85 


-1.00 


164Dy 


-0.93(9) 


0.347 


7(7) 

V / 


-0.77 


-1.00 


i54Gd 


-1.00(3) 


0.310 


1(5) 


-0.77 


-0.84 


i52Sm 


-0.98(4) 


0.307 


4(4) 


-0.54 


-0.83 


1880s 


-0.63(5) 


0.185 


17(3) 


-0.65 


-0.76 


190Os 


-0.35(9) 


0.177 


23(3) 


-0.53 


-0.75 


1920s 


-0.3(1) 


0.167 


25(2) 


-0.48 


-0.71 


194pt 


0.53(4) 


0.143 


41(1) 


0.08 


0.14 


196pt 


0.59(7) 


0.129 


42(2) 


0.00 


0.02 


106pd 


-0.4(1) 


0.230 


22(3) 


-0.41 


-0.55 


108pd 


-0.6(4) 


0.242 


17(16) 


-0.23 


-0.30 


112Cd 


-0.4(1) 


0.181 


22(2) 


-0.57 


-0.83 


ii^Cd 


-0.4(1) 


0.184 


23(3) 


-0.36 


-0.59 
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FIG. 2: Top panel: R^bm calculated over the whole IBM parameter space confirms a good ful- 
fillment of the second order approximation to K^. Middle panel: The two transitional legs for 
fixed values of x = — \/7/2 (left) and C = 1 (right). The approximation misses slightly the phase- 
/shape-transitional parameter region. Bottom panel: Effective 7-deformations calculated from 
and K^^^^' calculated over the whole parameter range. All calculations are for N = 10 bosons. 
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FIG. 3: R^B^M ■> analog to the top panel of Figure 2, but using only the first order approximation. 



Deviations from unity are much larger than in the second approximation. 



16 




17 



